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Abstract. Using a purely analytic approach to gaseous and dark matter halos, we study the cross-correlation 
between the Sunyaev-Zel'dovich (SZ) sky and galaxy survey under flat sky approximation, in an attempt to acquire 
the redshift information of the SZ map. Fhe problem can be greatly simplified when it is noticed that the signals 
of the SZ-galaxy correlation arise only from hot gas and galaxies inside the same massive halos (i.e. clusters), and 
field galaxies make almost no contribution to the cross-correlation. Under the assumption that both the hot gas 
and galaxies trace the common gravitational potential of dark halos, we calculate the expected cross SZ-galaxy 
power spectra for the WMAP/Planck SZ maps and the SDSS galaxy sample at small scales 100 < I < 1000. It 
turns out, however, that it is not presently feasible to measure such small angular cross power spectra because 
of the high noise levels at I > 400 with the WMAP/Planck experiments. Future SZ observations with better 
angular resolutions and sufficiently wide sky coverages will be needed if this technique is applied for the statistical 
measurement of redshift distribution of the SZ sources. 
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1. Introduction 

Most of the baryons in the local universe exist in the form 
of warm-hot intergalactic medium with temperature of 
T ~ 10 5 — 10 7 K as a result of gravitationally driven shocks 
and adiabatical compression when they fall into large-scale 
structures and collapsed dark matter halos (e.g. Cen & 
Ostriker 1999; Dave et al. 2001). In the latter case, the 
very hot baryons gravitationally bound in massive halos 
such as groups and clusters usually manifest themselves by 
strong diffuse X-ray sources in terms of bremmstrahlung 
emission, which are directly detectable with current X-ray 
instruments. Moreover, the energetic electrons of the hot 
gas also interact the passing cosmic microwave background 
(CMB) photons through the so-called Sunyaev-Zel'dovich 
(SZ) effect, giving rise to a subtle change in the CMB spec- 
trum. Current SZ measurements of known clusters at high 
signal-to-noise are now routine (e.g. Carlstrom et al. 2002) 
and recent detection of the excess power relative to pri- 
mordial CMB anisotropy at arcminute scales has been suc- 
cessfully attributed to the statistical signals of the thermal 
SZ effect (Mason et al. 2003; Bond et al. 2003; Komatsu 
& Seljak 2002). 

However, both X-ray and SZ measurements contain 
no information about the redshift distributions of the 
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hot baryons and their host groups/clusters. Spectroscopic 
follow-up observations should be made to complement our 
knowledge of the location and evolution of the hot baryons 
and their host halos. It has been realized that a more 
practical and powerful approach to extracting the dis- 
tance information from the X-ray and SZ maps is per- 
haps to cross-correlate the X-ray and SZ maps with the 
existing galaxy catalogs or ongoing deep galaxy surveys 
(e.g. Seljak, Burwell & Pen 2000; Zhang & Pen 2001;2003; 
Zhang, Pen & Wang 2002). In particular, the problem 
can be greatly simplified if we notice that the signals of 
the cross-correlation between the SZ (or X-ray) map and 
galaxy survey arise primarily from cluster galaxies, and 
field galaxies make almost no contribution to the SZ(or 
X-ray)-galaxy correlation. This permits a straightforward 
calculation of the power spectrum of the SZ(or X-ray)- 
galaxy correlation from an analytic model of dark halo 
abundance, i.e., the Press & Schechter (1974; PS) formal- 
ism, along with a reasonable prescription of galaxy and 
gas distributions inside a given halo of mass M at red- 
shift z. For the latter, one can adopt either the simplest 
self-similar NFW-like profile suggested by numerical sim- 
ulations (Navarro, Frenk & White 1995; NFW), or the 
empirically motivated density profiles such as the King 
model or /3-model, in combination with the halo occupa- 
tion distribution (see Cooray & Sheth 2003 for a recent 
review). In this paper, we explore the power spectrum of 
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the SZ-galaxy cross-correlation based purely on the halo 
approach. Similar technique has been recently applied by 
Zhang & Pen (2003) to the study of the cross-correlation 
of the soft X-ray background and galaxy survey and by Wu 
& Xue (2003) to the study of the auto-correlation of the 
soft X-ray background. In their early work, Zhang & Pen 
(2001) also investigated the SZ-galaxy cross-correlation 
using the so-called continuum field model, in which the 
gas distribution in dark halos is given by a convolution 
of the dark matter distribution with a Gaussian window 
function, while the matter fluctuation in the highly non- 
linear regime (e.g. clusters) is calculated in terms of the 
hyper-extended perturbation theory. 

For the purpose of actual applications, we will calcu- 
late the expected cross-correlation between the SZ maps 
observed by WMAP/Planck and the galaxy survey by 
the Sloan Digital Sky Survey (SDSS), and assess the 
feasibility of extracting the redshift information of the 
WMAP/Planck SZ maps using this statistical approach. 
Throughout this paper we adopt a flat cosmological model 
(ACDM) with the best fit parameters determined by 
WMAP: Q A = 0.73, fl M = 0.27, h = 0.71, Q b h 2 = 0.0224, 
o- 8 = 0.84 and n s = 0.93 

2. SZ-galaxy cross-correlation 

2.1. Properties of gaseous and dark matter halos 

The thermal SZ effect along the direction 9 due to the hot 
gas inside a halo can be evaluated following 
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where N e and T arc the number density and tempera- 
ture of electrons, x = /ipf/feTcMB is the dimcnsionlcss 
frequency, Tcmb is the temperature of the present CMB, 
and the integral is performed along the line of sight. We 
specify the gas temperature in terms of virial theorem 
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where A c is the overdensity of dark matter with respect to 
the critical value Pcrit, which is calculated using a fitting 
formula given by Bryan & Norman (1998) for fx = 0.8, 
and E 2 {z) = fi M (l + zf + Q A . 

For simplicity we assume that both gas and galaxies 
in massive halos trace the dark matter distribution such 
that 
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in which /(, is the universal baryon fraction, /j, e — 1.13 is 
the mean electron weight, AA ga i(r) is the number density 



of galaxies, and A is a proportionality constant. Indeed, 
this so-called self-similar model provides a good approxi- 
mation for the radial profiles of gas and galaxies in dark 
halos if our focus is not on the detailed physical pro- 
cesses of the gas and galaxies. Unlike thermal X-ray emis- 
sion which is governed by the central gas content of mas- 
sive dark halos such as groups and clusters of galaxies, 
the thermal SZ effect is generated by the hot gas dis- 
tributed over whole halo regions. Beyond the core regions 
the radial profiles of gas and dark matter look roughly 
very similar in shape. In recent years some physically mo- 
tivated models have been proposed for the density and 
temperature of the hot gas in massive halos (e.g. Babul 
et al. 2002; Voit et al 2002; Wu & Xue 2002; etc.). One 
can easily extend our formalism to inclusion of these gas 
profiles. Here we adopt the universal density profile as 
suggested by NFW to describe the dark matter distribu- 
tion inside halos, Ap DM (r) = S ch p clit /[{r/r s )(l + r/r s ) 2 ], 
and fix the free parameter using the empirical fitting for- 
mula found by numerical simulations (Bullock et al. 2001): 
c = [9/(1 + z)](M/2.1 x 10 13 M Q )-°- 13 , where c = r vir /r s 
is the so-called concentration parameter. 

In order to determine the constant A in equation (6), 
we need to know the mean galaxy number for a given 
halo of mass M, i.e. the halo occupation distribution, 
Afgai(M). We employ the best fit analytic formulae of 
Sheth & Diaferio (2001) for the spiral [JV ga i iS (M)] and el- 
liptical [N ga \^{M)\ galaxies based on the GIF simulations 
(Kauffmannn et al. 1999): 

AW(M) = (M/M s )^+O.5 e - 4[log(M/loll - 75M )l2 ; 

AWe(M) = (M/M E )^ e - (2xl ° llM /M)2 ; 

iVgai(M) = JV ga i,s(M)+JV ga i, E (M) ) (7) 

where M s = 7 x lO 13 /i- 1 M , a s = 0.9, M E = 3 x 
1O 12 /i _1 M0, and cue = 0.75. The mean number density 
of galaxies is thus 
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dMdV 

where d 2 N/dMdV is the comoving number density of dark 
halos (see below). Given the halo occupation distribution, 
we are able to normalize the number density profile of 
galaxies within a halo of mass M following 



iVgai(M) = J AN &al {r)dV. 



(9) 



The integration above is made over whole cluster region 
out to virial radius. This gives the proportionality con- 
stant A in equation (6) 

47r5 ch pcritr s 3 [ln(l + C )- C /(l + c)]' k ' 

Now, one can easily get the surface number density of 
galaxies by projecting AA ga i(r) along the line-of-sight 6: 



n ga i(M, z 



6) = D\(z) J AN gal (r)d X , 
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in which Da(z) is the angular diameter distance to the 
halo. 
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2.2. Power spectra 

The angular cross power spectrum of the SZ-galaxy cor- 
relation can be separated into the Poisson term Cf and 
the clustering term Cf: 



Cf = g(x) J 
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where Dq is the comoving angular diameter distance to 
the halo of mass M at z, and yi and ni are the Fourier 
transforms of the Compton y-parameter and the surface 
overdensity of galaxies [n ga i(M, z, 0) - n gal (z )]/n gal (z ), 
respectively, b(M, z) is the bias parameter, for which we 
use the analytic prescription of Mo & White (1996), and 
the mean surface number density of galaxies n ga i(zo) is 
given by 



^gai(^o) 



N gal (z)(l + z) 2 B(z,z a )D 2 A (z) ^dz, 
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where and also in equations (12) and (13) B(z, z ) de- 
notes the redshift distribution function of galaxies. Here 
we use the following approximate expression (e.g. Baugh 
& Efstathiou 1993; Dodelson et al. 2002): 



B(z,z ) 
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in which z is the median redshift of galaxy distribution 
in galaxy survey. 

Finally, we adopt the PS mass function for the comov- 
ing number density of dark halos 
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in which S c is the linear over-density of perturbations that 
collapsed and virialized at redshift z, and a{M) is the 
linear theory variance of the mass density fluctuation in 
sphere of mass M: 
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and W(x) = 3(sinx — xcosi)/i 3 is the Fourier repre- 
sentation of the window function. The power spectrum, 
P(k) cx k ns T 2 (k), is normalized by the rms fluctuation on 
an 8 h^ 1 Mpc scale, ag, and we take the transfer function 
T(k) from an adiabatic CDM model given by Bardeen et 
al. (1986). If we replace ni by g{x)yi and set B(z, zq) = 1 in 



equations (12) and (13) we can get the thermal SZ power 
spectrum, C ; sz ~ sz . Instead, if g{x)yi is replaced by n;, we 
will obtain the galaxy-galaxy correlation, (7f a _sa . Note 
that for the later case the square of \m\ in equation (12) 
should be \m\ if (iV ga i(iV ga i - 1)) < 1 (Seljak 2000). 



3. Application: prediction for the cross-correlation 
between the MAP/Planck SZ maps and SDSS 
galaxy survey 

First, we demonstrate in Figure 1 the expected power 
spectra of auto-correlations of the SZ sky and galaxy 
survey and their cross-correlation at frequency v = 30 
GHz and for all galaxies (B = 1). Note that we have 
dropped the negative sign in g(x). Essentially, the SZ- 
SZ power spectrum is dominated by Poisson distribution 
over all multipole range from I = 100 (e.g. Komatsu 
& Kitayama 1999), while clustering of the sources gov- 
erns the galaxy-galaxy correlation at angular scales up 
to I w 5000. Overall, both of the SZ-SZ and SZ-galaxy 
power spectra exhibit a decline at sufficiently large I, re- 
gardless of the fact that the galaxy-galaxy correlation still 
has much power at high I beyond I = 10 4 . The latter 
can be attributed to the pairs of galaxies residing within 
the same, smaller halos at higher redshifts which arc in- 
cluded in the estimate of galaxy-galaxy correlation be- 
cause of the choice, B = 1. Extrapolation of the halo 
occupation number to high redshifts may lead to an over- 
estimate of galaxies in small halos. Note that unlike the 
Compton y-parameter, the galaxy surface number den- 
sity n ga i(M, z, 0) is proportional to D 2 A {z) though both of 
them have the same spatial dependence on the underlying 
dark matter distribution of clusters. The strong power of 
the galaxy-galaxy correlation at high I can be significantly 
suppressed when we work with a realistic redshift distri- 
bution function B(z,z ) based on galaxy surveys such as 
SDSS (see below). 

Next, following the investigation of Peiris & Spergel 
(2000), we carry out our numerical predictions for the 
galaxy sample of SDSS and the SZ sky (to be) revealed by 
WMAP/Planck surveyor. We study the feasibility of ex- 
tracting the redshift information of the SZ sky from their 
angular cross power spectra. Toward this end, we select 
the SDSS galaxies in different magnitude ranges charac- 
terized by different values of median redshift zq in the 
distribution function B{z,zq) (Dodelson et al. 2002). We 
assume that all the galaxies have photometrically deter- 
mined redshifts to optimize the utilization of the SDSS 
data. This yields a total of 1.2 x 10 s galaxies over a quar- 
ter of the sky in the magnitude bin 20 < r* < 21. We 

have also compared the angular power spectra of galaxies 
^gal-gal predicted by 

our halo model with those deter- 
mined by Tcgmark et al. (2002) from early SDSS data 
on scales I < 600 in four different magnitude ranges from 
r* = 18 to 22, and found that the agreement is good when 
the bias factors suggested by Tegmark et al. (2002) are 
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Fig. 1. The angular power spectra of the SZ-SZ (top 
panel), galaxy-galaxy (middle panel) and SZ-galaxy (lower 
panel) correlations at frequency v = 30 GHz and for all 
galaxies. Poisson and clustering components are explicitly 
displayed. 

used in different magnitude bins [see Fig. 2 of Tegmark et 
al. (2002)]. 

The working frequencies of both WMAP and Planck 
are fixed at v « 30 GHz. Whether or not one is able 
to effectively measure the SZ-SZ, galaxy-galaxy and SZ- 
galaxy power spectra depends critically on the noise levels. 
Following the standard treatment, we estimate the errors 
on the SZ-SZ power spectrum, C ; , and galaxy-galaxy 
power spectrum, C, ga _ga , through 
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sky 



where f^ y and are the MAP/Planck and SDSS sky 
coverages, respectively, and Al is the bin width that will 
be fixed to be 1/4 in our estimation (see Zhang et al. 2002). 
C* ot denotes all contributions to the measured power spec- 
tra that are essentially composed of true signals and detec- 
tor noise, C\ ot = Cf uo + C, noiso , provided that foreground 
contamination is properly removed. For the SDSS survey, 
the shot noise per mode is simply 



cr 



i 

n ga \ 



(20) 



For the SZ sky observed by MAP/Planck, the instrumen- 
tal noise per mode can be modeled by 



= w~ 1 {v)e l 



(21) 



in which we have assumed that the experimental beam 
is Gaussian with width 6 2 {v) [FWHM= (81n2) 1 / 2 0(^)], 
and w~ x / 2 denotes the noise level in CMB temperature 
fluctuation per pixel. For our working frequency of v « 30 
GHz, we choose these parameters in terms of the tabulated 
values given by Cooray & Hu (2000). In a similar way, we 
can estimate the error for the SZ-galaxy power spectrum 
from 



AC, 
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We show in Figure 2 the expected power spectra of 
the auto-correlations of the WMAP /Planck SZ maps and 
the SDSS galaxy survey and their cross-correlations for 
z Q = 0.33, corresponding to 20 < r* < 21. While in prin- 
ciple the galaxy power spectrum can be well constructed 
with the SDSS survey over a broad range of multipoles, the 
measurements of the SZ power spectra beyond I ~ a few 
hundreds are hampered by the poor angular resolutions 
of the WMAP/Planck experiments, although the result 
of Planck is slightly better [see Bennett et al. (2003) for 
the first year WMAP result] . Recall that on the observed 
CMB spectrum the SZ signals appear to be dominant only 
at > 2000 unless the primary CMB can be accurately re- 
moved. Consequently, it may become hopeless to measure 
the cross-correlation between the WMAP /Planck SZ sky 
and the SDSS galaxy survey. Yet, another possibility is 
to search for the cross-correlation signals at very large 
angular scales generated by superclusters (Hernandez- 
Monteagudo & Rubino-Martin 2003) or large-scale matter 
inhomogeneities at low redshifts (e.g. z < 0.1) (Afshordi, 
Loh & Strauss 2004). Indeed, the noise appears to be rel- 
atively smaller at low redshifts and large angular scales 
(see Figure 3), which facilitates the detection of the cross- 
correlation of the WMAP/Planck SZ maps with galaxy 
surveys such as SDSS and 2MASS, although in this case 
our primary goal of extracting rcdshift information of the 
SZ map from the cross-correlation signals becomes less 
interesting. 

Taking the Planck experiment as an example, we have 
shown in Figure 3 the SZ-galaxy cross-correlations for a 
set of median redshifts of the SDSS galaxies. Despite the 
large errors at high / beyond ~ 400, it demonstrates the 
actual procedures that one needs to perform in order to 
extract the redshift information of the SZ map from these 
cross-correlations. 

Finally, we calculate the correlation coefficient between 
the SZ map and galaxy survey following 



r(0 = 
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Fig. 2. Expected angular power spectra of the SZ-SZ (top 
panel), the SDSS galaxy-galaxy (middle panel) and the 
SZ - SDSS galaxy (lower panel) correlations at frequency 
v = 30 GHz and for z = 0.33. Bin width is chosen as 
Al = I /A. Left and right panels correspond to the WMAP 
and Planck results, respectively. 




Fig. 3. Dependence of the Planck SZ-galaxy power spectra 
on redshifts of the SDSS galaxies. 



Recall that r COII (l) acts as a quantitative indicator of the 
significance of the SZ-galaxy correlation. The associated 
error can be estimated by 



(Ar) 2 = r c 2 orr (0 
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We display in Figure 4 the resulting correlation coefficients 
in the same redshift ranges as in Figure 3 for the Planck 
SZ-SDSS galaxy correlation. Similarly to the power spec- 
tra in Fig. 3, instrumental noise becomes to be dominant 
at both lower and higher I for the correlation coefficients. 
Peak location of r corr increases with the median redshift 
of SDSS galaxy distribution because inclusion of more dis- 
tant galaxies in the estimate of the SZ-galaxy correla- 
tion raises the power at smaller angular scales. Overall, 
the correlation efficient has a value of r corr w 0.2 at 
100 < I < 1000, indicating that the two phenomena are 
only moderately correlated with each other. This arises 
from the fact that only can galaxies in clusters contribute 
to the SZ-galaxy correlation and field galaxies selected in 
SDSS sample have almost no effect on the cross corre- 
lation between the SZ map and galaxy survey. Zhang & 
Pen (2001) made a similar analysis to maximize r corr by 
properly picking the redshift distribution function of the 
SDSS galaxies, which can increase the correlation efficient 



to r" 



0.7 at a smaller scale of 10 3 < I < 10 4 . 



4. Discussion and conclusions 

Cross-correlation between the SZ map and galaxy survey 
will reveal, in a statistical manner, important information 
on the redshift distribution of hot baryons in the universe, 
which requires no spectroscopic follow-up observation of 
individual sources. Recall that the major advantage of 
measurement of the SZ power spectrum over the SZ clus- 
ter survey is that one can acquire weak SZ signals at high 
statistical significance level without the need for resolving 
individual clusters. Indeed, time-consuming spectroscopic 
follow-up observations of the SZ sources may eventually 
throw a shadow on the effective applications of the SZ 
power spectrum if redshift distribution is ultimately con- 
cerned. Cross-correlation between the SZ map and galaxy 
survey provides a simple approach which is just suited for 
the problem. In particular, theoretical prediction of the 
SZ-galaxy power spectrum can be greatly simplified if one 
notices that the signals of the SZ-galaxy cross-correlation 
are primarily produced by the hot gas and galaxies within 
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Fig. 4. The correlation coefficients for the Planck SZ map 
and SDSS galaxy survey in different redshift ranges. Bin 
width is taken to be Al = 1/4. 

the same clusters. Namely, field galaxies make almost no 
contribution to the SZ-galaxy cross-correlation. 

Assuming that both the intracluster gas and cluster 
galaxies follow the same dark matter distribution with a 
functional form of the NFW-like profile and the abun- 
dance of dark halos is described by the PS mass func- 
tion, we have predicted the SZ-galaxy angular power spec- 
trum. As it is expected, the power spectrum indeed shows 
a moderately strong correlation indicated by an overall 
correlation coefficient of r corr 0.2 at 100 < 1 < 1000. 
However, applying our algorithm to the SDSS galaxy sur- 
vey and MAP/Planck SZ maps yields an unpleasant re- 
sult: It is unlikely that one can acquire meaningful infor- 
mation about the redshift distribution of the MAP /Planck 
SZ maps because the SZ signals are significantly below the 
noise levels in our interested range of multipolcs, I > 1000, 
although the SDSS galaxy sample should be just suited for 
such a purpose. One possibility is to work at large angu- 
lar scales / < 100, which is related to the SZ effect pro- 
duced by nearby superclusters and clustering of galaxies 
(Hernandez-Monteagudo & Rubino-Martin 2003; Afshordi 
et al. 2004). Another application is to utilize the CMB 
data at very small angular scales beyond I = 2000 to be 
obtained by many ongoing/upcoming experiments such 
as ACBAR, AMiBA, BIMA, CBI, etc. In this case, one 
needs instead to deal with the small-sky coverage problem. 
Indeed, a similar work should be done on the feasibility 
of cross-correlating these small-sky coverage CMB data 
with existing galaxy catalogs. We conclude that the cross- 



correlation between the SZ map and galaxy survey can 
in principle yield valuable information about the redshift 
distribution of the host baryons in the universe. However, 
the actual applications of this technique to real observa- 
tions may not be possible until high angular resolution and 
sensitivity SZ power spectrum over a wide sky coverage is 
achieved. 
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